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Abstract.  The  currently  adopted  strategy  of  EOP  estimation  from 
VLBI  is  to  estimate  six  parameters:  UT1,  UT1  rate,  pole  positions  and 
nutation  offsets  for  each  24-hour  session  independently.  Then  the  result¬ 
ing  time  series  of  raw  Earth  orientation  parameters  are  filtered  and  a 
regression  analysis  is  performed  in  order  to  obtain  nutation  coefficients, 
polhode  of  the  pole  and  other  physical  parameters.  Thus,  the  latter  pa¬ 
rameters  are  obtained  indirectly  in  two  stages.  An  alternative  approach 
of  direct  estimation  of  the  final  Earth  orientation  parameters  is  presented. 
Pole  coordinates  and  UT1  are  considered  as  a  sum  of  three  components: 
the  low-period  component  which  is  modeled  by  a  cubic  spline,  the  har¬ 
monic  component  which  includes  forced  nutation,  precession  and  sub¬ 
daily  variations  of  EOP,  and  the  stochastic  component  which  is  modeled 
by  a  linear  spline  with  segment  length  1-2  hours.  All  parameters  are 
obtained  in  a  single  LSQ  solution  using  all  available  data. 


1.  Introduction 

Rotation  of  the  Earth  is  a  complicated  process.  We  don’t  observe  it  directly,  but 
estimate  parameters  of  the  model  of  indirect  observations.  The  choice  of  the  set 
of  parameters  is  inspired  by  both  theory  and  our  ability  to  estimate  parameters 
of  the  model.  Earth  rotation  in  the  reference  system  tied  to  the  rotating  Earth 
can  be  represented  by  a  superposition  of  three  components:  1)  low-frequency 
variations  (period  >  2d),  2)  quasi-diurnal  and  quasi- semidiurnal  variations,  and 
3)  a  stochastic  component. 

The  traditional  approach  of  EOP  estimation  from  VLBI  is  first  to  estimate 
six  parameters:  pole  coordinates,  UT1  and  its  rate,  nutation  offset  angles  for 
each  24-hour  session  independently,  then  to  filter  the  raw  time  series  of  the 
estimates,  and  then  to  perform  statistical  analysis  of  the  resulting  time  series. 
High-frequency  Earth  orientation  variations  cannot  be  determined  in  this  way, 
and  they  are  obtained  in  a  different  solution. 

The  traditional  approach  has  an  internal  contradiction.  If  we  assume  that 
nutation  is  a  purely  harmonic  process,  with  known  frequencies,  which  corre¬ 
sponds  to  a  quasi-diurnal  polar  motion  then  we  need  to  estimate  only  ampli¬ 
tudes  and  phases  of  the  harmonic  components,  but  not  to  represent  nutation 
angles  as  a  time  series.  If  we  assume  that  there  are  non-harmonic  components 
in  nutation,  e.g.  free  core  nutation,  with  a  broad  spectrum,  then  the  traditional 
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procedure  cannot  strictly  separate  them  from  pole  coordinates.  The  deficiencies 
of  the  traditional  approach  are: 

•  The  components  of  the  retrograde  quasi-diurnal  polar  motion  (nutation) 
with  frequencies  rather  far  away  from  the  diurnal  frequency  are  not  per¬ 
fectly  separated  from  the  long-period  components  of  polar  motion.  As  a 
result  the  estimates  of  their  amplitudes  are  corrupted. 

•  The  correlations  between  the  estimates  of  daily  nutation  offsets  are  ne¬ 
glected,  and  series  are  usually  treated  as  statistically  independent. 

•  Nutation,  precession,  high-frequency  EOP  variations  are  treated  differ¬ 
ently. 


2.  Direct  approach 

An  alternative  way  is  to  estimate  all  components  of  EOP  in  a  single  global 
solution.  Since  the  length  of  the  data  set  of  regular  VLBI  observations  started 
in  January  1984  is  increasing,  it  is  already  feasible  to  separate  aliases  related  to 
the  18.6-year  period.  Not  only  the  main  components  of  nutation  but  almost  all 
terms  can  be  determined  from  the  observations.  The  following  parameterization 
is  proposed: 

•  X p,  Yp ,  UT 1  and  their  rates  for  each  session.  They  model  slowly  varying 
components. 

•  Set  of  sin-  and  cos-  terms  for  Xp,  Yp  and  UT  1  at  frequencies  within  the 
diurnal  and  semidiurnal  bands  (with  respect  to  the  rotating  Earth).  They 
include  precession,  CEP  offsets,  nutations,  tidal  diurnal  and  semidiurnal 
EOP  variations,  harmonic  variations  at  ad  hoc  frequencies.  The  frequencies 
are  taken  from  theory. 

•  Parameters  of  a  linear  spline  for  Xp,  Yp  and  UT  1  with  segment  length  1.5 
hours.  They  model  the  stochastic  component  in  the  Earth’s  rotation. 

•  Other  parameters. 

Such  a  choice  of  parameterization  should  adequately  model  all  constituents 
of  Earth  rotation.  However,  these  parameters  create  some  difficulties: 

•  The  theoretical  spectrum  of  EOP  variations  has  closely-separated  compo¬ 
nents  (multiplets)  which  cannot  be  resolved  using  the  given  data  set. 

•  There  are  synonymous  constituents,  e.g.  nutation  and  diurnal  high-frequency 
variations  in  polar  motion. 

•  Three  constituents  are  overlapping  and  are  not  independent: 

—  there  is  a  linear  dependence  between  a  session-average  parameter  and 
an  average  value  of  linear  spline; 
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—  there  is  a  linear  dependence  between  a  session-rate  parameter  and  a 
rate  of  change  of  the  linear  spline; 

—  there  is  a  quasi-linear  dependence  between  a  harmonic  parameter  and 
coefficients  of  the  linear  spline. 

Fortunately,  we  can  overcome  these  problems.  The  lowest  critical  frequency 
is  assumed.  The  harmonic  variations  with  the  frequency  lower  than  the  critical 
frequency  are  considered  as  a  secular  change  and  they  are  represented  by  a  linear 
model.  The  harmonic  constituents  with  frequencies  which  differ  by  less  than  the 
critical  frequency  are  grouped  into  multiplets.  It  is  assumed  that  the  frequencies 
of  the  components  of  a  multiplet  and  the  ratios  of  the  amplitudes  are  known 
from  theory.  The  complex  admittance  of  the  multiplet  is  estimated  instead  of  the 
amplitude  and  the  phase  of  a  singlet.  Strong  de- correlation  constraints  on  the 
coefficients  of  the  linear  spline  for  each  session  are  imposed:  a)  the  average  is  zero 
(de-shifting),  b)  the  average  rate  is  zero  (de-trending),  c)  the  amplitudes  of  sin- 
and  cos-  components  at  the  specified  frequencies  are  zero  (de-harmonization). 


3.  Solution 

The  approach  outlined  above  has  been  implemented  in  the  CALC/SOLVE  pro¬ 
gram  package.  I  used  all  available  VLBI  observations  for  1979.7-2000.0,  2.8 -106 
observations,  2961  sessions. 

The  estimated  Earth  orientation  parameters  were: 

•  736  harmonic  parameters  including  precession  rate  and  obliquity  rate,  ce¬ 
lestial  pole  offsets,  148  nutation  frequencies,  34  additional  tidal  frequencies, 
3  ad  hoc  frequencies.  The  lowest  frequency  was  1  •  10-8  rad/sec  (period 
of  18.6  years) 

•  2559  daily  pole  coordinates  and  UT1  as  well  as  their  rates. 

•  44150  parameters  of  linear  splines  for  pole  coordinates  and  UT1.  The 
duration  of  a  span  between  nodes  of  the  spline  was  90  minutes. 

•  Other  parameters:  station  positions  and  velocities,  source  coordinates, 
troposphere  path  delay,  troposphere  gradients,  clock  parameters. 

Weighted  root  mean  square  of  the  postht  residuals  is  24.6  psec.  Formal 
uncertainties  of  the  estimates  of  harmonic  variations  are  20-300  prad.  The 
largest  correlation  of  0.91  is  between  the  18.6-year  nutation  term  and  precession. 
Only  89  correlations  out  of  270  480  are  greater  than  0.7. 

Comparison  with  the  empirical  nutation  coefficients  obtained  by  Souchay 
et  al.  (1995)  and  with  semi-empirical  coefficients  IERS96  (1996)  showed  the 
differences  are  small.  They  don’t  exceed  1-2  nrad,  but  they  are  not  negligible. 
Some  terms  are  presented  in  Table  1.  Plots  of  the  estimates  of  the  stochastic 
component  in  Earth  rotation  are  presented  in  Figure  3.. 
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Table  1.  Differences  in  some  terms  with  respect  to  the  IAU  1976 
precession,  IAU  1980  nutation  models. 


Term 

D 

irect  solution 

Souchay  (1995) 

IERS96 

Units 

Precession  rate 

-2.95  ±  0.02 

-3.00 

-2.96 

mas/yr 

Obliquity  rate 

—0.27  ±  0.009 

-0.26  ±0.01 

-0.23 

mas/yr 

AT  (in)  Per  18.“6 

-6.71  ±  0.07 

-7.63  ±0.13 

-6.67 

mas 

AT  (out)  Per  18.“6 

3.76  ±  0.11 

2.91  ±0.07 

3.65 

mas 

Ae  (in)  Per  18.“6 

2.64  ±  0.05 

2.86  ±0.03 

2.86 

mas 

Ae  (out)  Per  18.“6 

1.37  ±  0.03 

1.49  ±0.05 

1.55 

mas 

AT  (in)  Per  182.d6 

1.63  ±  0.02 

1.64  ±0.04 

1.69 

mas 

AT  (out)  Per  182.d6 

—  1.19  ±  0.02 

-1.42  ±0.04 

-1.40 

mas 

Ae  (in)  Per  182.d6 

-0.60  ±  0.01 

-0.55  ±0.01 

-0.54 

mas 

Ae  (out)  Per  182.d6 

—0.24  ±  0.01 

-0.46  ±0.01 

-0.46 

mas 

AT  (in)  Per  13.d66 

-0.27  ±  0.02 

-0.31  ±0.04 

-0.32 

mas 

AT  (out)  Per  13.d66 

0.20  ±  0.02 

0.34  ±0.04 

0.27 

mas 

Ae  (in)  Per  13.d66 

0.131  ±  0.007 

0.15  ±0.01 

0.164 

mas 

Ae  (out)  Per  13.d66 

0.127  ±  0.007 

0.17  ±0.01 

0.136 

mas 

4.  Conclusion 

It  is  shown  that  an  empirical  expansion  of  harmonic  variations  in  Earth  rotation 
can  be  obtained  in  one  step. 

•  The  advantages  of  the  proposed  approach  are 

—  A  uniform  treatment  of  harmonic  variations  in  the  Earth  rotation  at 
diurnal  and  semidiurnal  frequencies; 

—  No  Aip,  Ae  time  series; 

—  An  ability  to  estimate  correctly  short-period  nutation  amplitudes. 

•  The  disadvantage  is  high  computation  cost. 
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Figure  1.  Residual  variations  in  the  Earth  rotation  for  a  routine 
VLB  I  experiment 


